First principles elastic constants and electronic structure of Q:-Pt2Si and PtSi 
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We have carried out a first principles study of the elastic properties and electronic structure 
for two room-temperature stable Pt silicide phases, tetragonal Q-Pt2Si and orthorhombic PtSi. We 
have calculated all of the equilibrium structural parameters for both phases: the a and c lattice 
constants for a-Pt2Si and the a, b, and c lattice constants and four internal structural parameters 
for PtSi. These results agree closely with experimental data. We have also calculated the zero- 
pressure elastic constants, confirming prior results for pure Pt and Si and predicting values for the 
six (nine) independent, non-zero elastic constants of a-Pt2Si (PtSi). These calculations include a full 
treatment of all relevant internal displacements induced by the elastic strains, including an explicit 
determination of the dimensionless internal displacement parameters for the three strains in a-Pt2Si 
for which they are non-zero. We have analyzed the trends in the calculated elastic constants, both 
within a given material as well as between the two silicides and the pure Pt and Si phases. The 
calculated electronic structure confirms that the two silicides are poor metals with a low density 
of states at the Fermi level, and consequently we expect that the Drude component of the optical 
absorption will be much smaller than in good metals such as pure Pt. This observation, combined 
with the topology found in the first principles spin-orbit split band structure, suggests that it may 
be important to include the interband contribution to the optical absorption, even in the infrared 
region. 

PACS 62.20.Dc, 71.20.Be, 71.15.Nc 



I. INTRODUCTION 



Metallic Pt silicide compounds are used to make recti- 
fying junctions on silicon substrates. The Schottky bar- 
rier, which determines the activation energy for trans- 
port of the charge carriers (holes), i|a|-220-240 meV for 
orthorhombic PtSi on p-type Si(001),EfH matching an im- 
portant atmospheric "transparency window" in the in- 
frared region. For this reason these materials are well 
suited to infrared detector applications. Orthorhombic 
IrSi has also been used in these applicatiops and has a 
Schottky barrier of 160 meV on Si(001),Q matching a 
lower portion of the transparency window. An impor- 
tant advantage of the silicides is that they are more com- 
patible with current silicon-based fabrication technology 
than infrared sensitive semiconductors such as InSb and 
Hg,Cdi_xTe. 

Despite their importance in infrared detector appli- 
cations, relatively little theoretical work has been done 
to investigate the fundamental electronic structuce and 
equilibrium properties of the Pt silicides. Bisi et al.u used 
the iterative extended Huckel method to calculate the 
angular-momentum-resolved density of states (DOS) for 
a number of near noble-metal silicides, including a-Pt2Si 
and PtSi. Yarmoshenko et aln compared x-ray photo- 
electron spectroscopy (XPS) and x-ray emission spec- 
troscopy (XES) measurements to the electronic struc- 
ture calculated using a linear mufhn-tin orbital (LMTO) 
method for a number of 4d and 5d silicides, including 
PtSi. In order to provide a more complete understand- 



ing of the properties of the Pt silicides we have carried 
out first principles electronic structure calculations for 
two room-temperature stable phases, tetragonal a-Pt2Si 
and orthorhombic PtSi. In addition to calculating all of 
the equilibrium structural parameters (Sec. ||) we have 
also obtained values for all of the zero-pressure elastic 
constants (Sec. HI) for both phases. We have investi- 



gated the electronic structure (Sec. IV) which is directly 
relevant to the infrared detector applications of these ma- 
terials. The purpose in all of these calculations has been 
to provide a fundamental understanding of the ground 
state properties of the Pt silicides. We summarize our 
results in Sec. 0. 



II. ATOMIC STRUCTURE 
A. Crystal structures 

Table || summarizes the characteristics of the equilib- 
rium crystal structures for each of the materials consid- 
ered here. The pure Pt and Si constituents of the Pt sili- 
cides both crystallize in a cubic structure under normal 
conditions, face.centered cubic (fee) for PtQ and cubic 
diamond for Si.Ll Both cubic structures are characterized 
by a single lattice constant a (Table ||) . 

The conventional unit cells of the two Pt silicides a- 
Pt2Si and PtSi are illustrated in Fig. |^. The structure of 
the room-temperature (T < 968 K) a-phase of Pt2Si is 
body-centered tetragonal (bet) and resembles a distorted 
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CaF2 structureBQ A central Si atom is surrounded by 
eight Pt atoms, which are located in the corners of a rect- 
angular cell elongated along the c-axis. There are two 
symmetry-equivalent Pt and one Si atom in the prim- 
itive cell. The unit cell is characterized by two lattice 
constants a and c (Table ||). 

PtSi has a primitive orthorhombic structure [see 
Fig. |l|(b)] with four symmetry-equivalent Pt and four 
symmetry-equival£iiL_Si atoms per primitive cell in an 
MnP-type lattice.EjO This structure is characterized by 
three lattice constants, denoted a, &, and c (Table ||). 
Half of the Pt and half of the Si atoms in the primitive 
cell are located in a (010) plane at with the rest of the 
atoms in a (010) plane at |6. The in-plane atomic coor- 
dinates are not completely specified by the space group 
symmetry and thus there are four free internal structural 
parameters upt, fpt, "Si: and wsi (the "1st position" col- 
umn in Table m. The experimental values are given in 



Table til. Each Si is surrounded by six Pt atoms at the 



corners of a distorted trigonal prism. The Pt atoms have 
six Si neighbors at the corners of a distorted octahedron, 
and four Pt neighbors, wJiich are positioned in four of 
the octahedral interstices.EJ The Pt atoms are arranged 
in zig-zag chains along the [100] direction. 



B. FPLMTO method 

The equilibrium structural parameters and zero- 
pressure elastic constants were calculated with a ful]j)p=. 
tential hnear muffin-tin orbital (FPLMTO) methodB'El 
which makes no shape approximation for the crystal 
potential. For mathematical convenience the crystal 
is divided into regions inside atomic spheres, where 
Schrodinger's equation is solved numerically, and an in- 
terstitial region. In our FPLMTO method the basis 
functions in the interstitial region are smoothed Han- 
kel functions.EJ This method does not require the use 
of empty spheres, even for open structures such as cubic- 
diamond-phase Si. The atoms were treated scalar rel- 
ativistically within the local density approximation, us- 
ing the-iexchange-correlation potential of Ceperley and 
Alder.E3 Spin-orbit interactions were not included. The 
choice of basis functions for the Pt and Si atoms was op- 
timized according to the procedure described in Ref. [l^. 
The parameters describing the basis are listed in Tablepv] 
and were used in the calculations for both silicides in ad- 
dition to pure Pt and pure Si.tS The Si 3s, 3p, 3d, and 
4/ as well as the Pt 6s, 6p, 5d, and 5/ were all included 
as valence orbitals. The Pt semi-core 5s and 5p were 
treated as core orbitals and we have not used the frozen 
overlapped core approximation (FOCA — see Ref. [T^ ). 

The equilibrium volume Vq and bulk modulus Bq of Pt 
and Si were determined by fitting the total energy calcu- 
lated at nine diffepjut, lattice constants to a Murnaghan 
equation of state.ElHia In the case of Pt (Si) we used a 
24x24x24 (12x12x12) cubic special k-point mesh which 



gave 6912 (864) points in the full Brillouin zone (BZ) and 
182 (28) points in the irreducible wedge. In addition, a 
real space mesh is used for calculating integrals of the po- 
tential over the interstitial region. We used a 16x16x16 
mesh for Pt and an 18x18x18 mesh for Si. With these 
choices the total energies were converged to better than 
1 /xRy per atom. 

In order to determine the equilibrium lattice constants 
of the silicides, the total energy hypersurfaces were mini- 
mized simultaneously with respect to all of the lattice pa- 
rameters. In the case of tetragonal Q!-Pt2Si, total-energy 
calculations were performed at nine different values of 
each of the two lattice constants a and c (a total of 81 
calculations) in the range 0.93 < a/aoxp, c/ccxp < 1.07. 
For orthorhombic PtSi we used seven different values of 
a, 6, and c, respectively (343 calculations), within the 
same relative ranges. For both Q;-Pt2Si and PtSi, 216 k- 
points were sampled in the full BZ which reduced to 28 
and 27 special k-points in the irreducible wedge, respec- 
tively. In the case of a-Pt2Si a 16 x 16 x 16 real space mesh 
was used for the interstitial integrals and a 24x16x24 
mesh was used for PtSi. With these choices the en- 
ergy per atom differed from the fully converged value 
(which was found for approximately 13 000 k-points in 
the full BZ) by 0.20 mRy for Q;-Pt2Si and 0.19 mRy for 
PtSi. The resulting total energy hypersurfaces, E{a, c) 
and E{a,b,c), were each fit to a third order polynomial 
in the lattice parameters. These polynomials were then 
minimized to yield the equilibrium parameters. The bulk 
modulus and equilibrium volume were obtained by fitting 
the total energy as a function of volume to a four-term 
Birch-Murnaghan equation of state,Ej 



(1) 



This fit was differentiated twice, yielding the bulk mod- 
al the theoretical equilibrium 
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Vo 



ulus Bo = V 
volume. 

The four free internal structural parameters of or- 
thorhombic PtSi, (wpt, vpt, usi and vsi in Table |), were 
determiaod sclf-consistently by calculating the ah initio 
forcesllSEj on tha,ions and, within the Born-Oppenheimer 
approximation,^ relaxing the position of each individual 
atom in the direction of the forces until the absolute val- 
ues of the forces were converged to less than 1.5 mRy/a.u. 
512 special k-points were used within the full BZ (cor- 
responding to 64 in the irreducible wedge). Initially the 
atomic positions-were relaxed starting from the experi- 
mental structurdiil and holding the three lattice constants 
fixed at their experimental values. Using these theo- 
retically determined internal parameters, the theoretical 
equilibrium lattice constants, oq, bo, and cq, as well as the 
bulk modulus Bo, were determined using the procedure 
described above. A second geometry relaxation was then 
carried out but now holding the lattice constants fixed at 
these theoretically determined equilibrium values. This 
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yielded a second set of internal structural parameters 
which we refer to as the self-consistent theoretical values. 
In principle this cycle could be repeated many times to 
obtain a set of lattice constants aad internal parameters 
which are truly "self-consistent. "E3 However, in practice 
we find that after the first cycle there are only small dif- 
ferences between the two sets of internal parameters (see 
Table ill) and so we regard them as being converged. We 
also note that the value of Bq obtained using fixed values 
of the internal structural parameters, as described here, 
is not strictly correct and that we relax the constraint 
of fixed internal parameters when we discuss the elastic 
constants (including Bq) in Sec. HIE. 



C. Equilibrium properties 



In order to test our method0 and, in particular, to 
test our choice of basis functions, we calculated the equi- 
librium lattice constant oq and bulk modulus Bq for Pt 
and Si as described in Sec. |II B| . The results are given in 
Table |l| and compared to experimental data. Since we 
will focus next on the elastic constants, we pay particular 
attention to Bq, which is essentially an elastic constant. 
We see in Table || that we obtain rather good agree- 
ment between our values and the experimental data. The 
self-consistent equilibrium lattice-constants gq and cq for 
Q!-Pt2Si are also listed in Table |l[ along with the theo- 
retical oo, bo, and cq lattice constants for PtSi. The bulk 
moduli Bq at the theoretical volumes for both Pt silicides 
are also given. 

In Tabic II] we list the internal structural parameters 
for PtSi. The values calculated using the experimental 
lattice constants are very close to those measured by 
Graeber et al\13 Comparing the atomic positions using 
the experimental versus the theoretical internal parame- 
ters, we find absolute shifts in the positions of less than 
0.028 a.u. The self-consistent internal parameters ob- 
tained using the theoretical lattice constants are gener- 
ally even closer to the experimental values. In this case 
the absolute atomic shifts relative to the experimental 
geometry are less than 0.023 a.u. 

The theoretical cohesive energies -Ecoh and heats of for- 
mation AH{ for the self-consistent equilibrium atomic 
geometries are compared to the experimental values in 
Table 0. As is typically the case, our local-density- 
functional based calculations overestimate the cohesive 
energy. However, the calculated heats of formation are 
much closer to the experiment. The experimental heats 
of formation are given for T = 298.15 K whereas the 
theoretical values correspond to K and do not include 
corrections for zero-point vibrations. We note that the 
heats of formation are very similar for the two silicides 
and that a plot of the theoretical AHf as a function of 
atomic percent Pt is concave up, as required for the sili- 
cides to both be thermodynamically stable. 



III. ELASTIC CONSTANTS 

A. Method of calculation 

The elastic constants determine the stiffness of a crys- 
tal against an externally applied strain. For small defor- 
mations we expect a quadratic dependence of the crystal 
energy E on the strain (Hooke's law). The elastic con- 
stants Cijki describe this quadratic behavior. Consider 
a displacement u(R) which takes every Bravais lattice 
point R of the undistorted lattice to a new position R' 
in the strained lattice, 



(2) 



where the index i corresponds to Cartesian coordinates. 
If we assume the applied strain is homogeneous (uniform 
throughout the crystal), we can rewrite Eq. |^ as 



with 



lij — Sij 



(3) 



Eor a homogeneous applied strain the displacement gradi- 
ents dui (R) /dRj are simply constants, independent of R. 
These displacement gradients define the nine components 
of a tensor. However, since the total energy E cannot 
change under rotations of the crystal as a whole, E can. 
only depend on the symmetric part of the deformation,Lj 
called the strain tensor e: 



dRi ^ 



OR, 



(4) 



Expanding the internal energy iJQ/,e) of the crystal with 
respect to the strain tensor givesEJ 



EiV,{e,nn})^EiV) + Vj2 
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where the stress tensor cr is defined by 
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(7i,j = — 
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dei 



(6) 



the second order adiabatic elastic constants are given by 
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Cijkl — y 



d^E{V,{e^^}) 
deijdeki 



(7) 



and V is the volume of the unstrained crystal. It is con- 
venient to use Voigt notation which takes advantage of 
the symmetries of the tensors: xx — > 1, yy ^ 2, zz — > 
3, yz — > 4, xz f—t 5 and xy 6. Using this notation 
Eq. d) becomesEl 



V ■ 



EiV, {ci}) = E{V) -f 1/^ (T,e, + ^Y. 



(8) 
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with the strain tensor given by 
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In order to calculate all M elastic constants of a crys- 
tal we applied M independent strains e*-^^ to the unit 
cell, using Eqs. ^ and to determine the atom po- 
sitions within the strained unit cell. In particular, we 
have M = 3 for both cubic Si and cubic Pt, M — 6 
for tetragonal Pt2Si, and M — 9 for orthorhombic PtSi. 
Each strain / = 1 , . . . , M was parameterized by a sin- 
gle variable 7 and we calculated the total energy E'^^\'^) 
for a number of small values of 7. For these small dis- 
tortions, i?'^' (7) was fit to a polynomial in 7 and then 
equated to the appropriate elastic constant expression 
E{V,{e^P in Eq. (§). From all of the fits we ob- 
tained a system of M linear equations for the elastic con- 
stants, which was solved for the Cij . Since we always take 
the undistorted crystal to be the zero-pressure theoreti- 
cal equilibrium structure, the applied stress cr is zero and 
so the second term of Eqs. (||) and (||) does not enter in 
the calculations described here. 

The parameterizations we used for the three indepen- 
dent strains in the cubic cases of Pt and Si are given in 
Table 0. Strain / = 1 is a volume-conserving stretch 
along the z-axis, the second strain is equivalent to sim- 
ple hydrostatic pressure, and strain / = 3 corresponds to 
a volume-conserving monoclinic shear about the z-axis. 
We carried out calculations for 9 values of 7 in the range 
of —0.01 to 0.01 for strains 1 and 2. However, for strain 
3 we calculated 9 points in the range from —0.04 to 0.04 
because the changes in the energy were rather small (a 
maximum of 0.1 mRy for 7 — 0.01), leading to larger er- 
ror estimates in the case of the smaller range. In order to 
calculate the six independent and non-vanishing elastic 
constant s of tetragonal a-Pt2Si we used the strains given 
in Table VII. E3 Orthorhombic PtSi has nine independent 
elastic constants and we chose the nine strains listed in 
Table VIII , For each of the silicide strains we carried out 
calculations for seven values of 7 in the range of —0.01 to 
0.01, except for strains 8 and 9 in the case of PtSi where 
only five values of 7 were considered (these monoclinic 
strains were particularly CPU intensive). Calculational 
errors in the elastic constants were determined from the 
least-squares fit to -E(7). All of our results were obtained 
from fits of the energy to third order in 7 because these 
yielded the smallest errors compared to polynomials of 
order two and four; the single exception was Pt C44, where 
minimum standard errors resulted from a fourth order fit. 

Calculations of the elastic constants require a very 
high degree of precision because the energy differences 
involved are of the order of 10 to 1000 /iRy. This cir- 
cumstance requires the use of a fine k-point mesh. With 
our choice of 23 328 special k-points in the full BZ for 
Pt, 864 for Si, and 5832 for a-Pt2Si and PtSi, the en- 



ergy per atom was converged to 1 /zRy or better in all 
cases. In order to minimize numerical uncertainties we 
used the same k-point mesh for all of the calculations in 
a given material. The differing symmetries of the vari- 
ous strains / resulted in differing numbers of irreducible 
k-points. We also checked that we obtained the same 
total energy for 7 = 0, regardless of strain / (and hence 
different symmetry and irreducible k-points). All of the 
calculations were carried out at the theoretical equilib- 
rium lattice constants listed in Table ||. Relaxation of 
the internal degrees of freedom was carried out in the 
case of all nine PtSi elastic constants. These relaxations 
are necessary because the atomic positions are not com- 
pletely fixed by the space group symmetry, even for the 
unstrained crystal, and consequently there exist free in- 
ternal parameters (see Table III) which must be rede- 
termined for any distortion of the crystal, including hy- 
drostatic pressure. Relaxations were also carried out in 
those cases where the strain-induced symmetry-reduction 
prompted it (C44 for Si and strains 1, 4, 6 for a-Pt2Si). 
For comparison we have calculated "frozen" elastic con- 
stants in these same cases, where the internal structural 
parameters where frozen at their zero-strain equilibrium 
values. 



B. Pt 



The three elastic constants for Pt are listed in Ta- 
ble IX. Pt is the only one of the metals considered in 
this work for which experirrietital data on elasticity is 
available. MacFarlane et aZ.Ea extrapolated the values 
to K, which makes them well suited for a comparison 
to our zero-temperature calculations. In the case of cn 
and C12 we find good agreement between our results and 
the experimental data (within 3-4%). The value of C44 
deviates by 14%, although the absolute error is approxi- 
mately 10 GPa for all three elastic constants. The error 
in C44 can be understood if we look closely at the band 
structure. Pt exhibits a wealth of van Hove singulari- 
ties directly at the Fermi energy, making it difhcult to 
integrate over the Fermi surface. A high density of k- 
points (23 328 in the full BZ) and a very small smearing 
width of 7 mRy in the higher-order smearing procedurecj 
are essential because the Fermi energy and hence the to- 
tal energy depend quite sensitively on these parameters. 
The calculated value of C44 was found to be more sensi- 
tive to the k-points than the other two elastic constants. 
Conversely, the silicidcs did not warrant such a special 
treatment and were calculated with a smearing width of 
25 mRy. It seems plausible that a more accurate treat- 
ment of the elastic properties of Pt may also require in- 
clusion of spin-orbit coupling. The bulk modulus calcu- 
lated from the theoretical values of the elastic constants 
[Bo = \{cii + 2ci2)] is 290.8 GPa. It agrees weU with 
both the experimental value of 288.4 GPa and the one 
extracted from the fit to a Murnaghan equation of state. 
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287.8 GPa (Sec. [nq ). 

The requirement of mechanical stabihty in a cubic 
crystal leads to the following restrictions on the elastic 
constantscj 

(cii - C12) > 0, cii > C44 > 0, (cii + 2ci2) > 0. 

(10) 



The Pt elastic constants in Table [X obey these stability 
conditions, including the fact that C12 must be smaller 
than cii. These conditions also lead to a restriction on 
the magnitude of Bq. Since Bq is a weighted average of 
cii and C12 and stability requires that C12 be smaller than 
cii, we are left with the result that Bq is required to be 
intermediate in value between cn and C12, 



C12 < Bo < cii 



C. Si 



(11) 



Because Si has a cubic structure, it has only three dis- 
tinct, non- vanishing elastic constants. These were deter- 
mined with the same strains as in the case of Pt (Ta- 
ble [v| ). Our results are close to experiment, as in- 
dicated in Table IX. The bulk moduli from the to- 



tal energy minimization and from the elastic constants 
[Bo = |(cii -f 2C12)] have the same value of 95.9 GPa, 
close to the one calculated from the experimental elastic 
constants, 97.0 GPa. The Si elastic constants in Table IX 
also obey the cubic stability conditions in Eq. (p^), mean- 
ing that C12 < Bo < cii. 

It is perhaps worth noting that the calculation of C44 
required a relaxation of the positions of the Si atoms 
within the distorted unit cell. The symmetry reduction 
by the monoclinic shear (e*^"^^ in Table VI) allowed the 
Si atoms to relax in the [001] direction. Without this 
relaxation, C44 would have been 108.6 GPa; this is to 
be compared with the relaxed value, of 79.9 GPa and 
the experimental value of 79.1 GPaEJ We have also ob- 
tained the dimensionless Kleinman internal displacement 
parameter ( which determines the magnitude of the in- 
ternal displacements along the [001] direction, 



Aw 



(/=3) 



(12) 



where a is the lattice constant and Cxij =^eg is the appro- 
priate element of the strain tensor [Eq. iM) and Table |Vlj . 

Fitting our calculated values of Aug^" to a quadratic 
function in 7 we find a value of C = 0.53 which agrees 
very well with the experimental value of 0.54.Ld 



D. a-PtaSi 



We have applied the six strains listed in Table VII in 



order to determine the clastic constants of tetragonal a- 
Pt2Si. The orthorhombic strains 1 and 4, and the mon- 
oclinic strain 6 all reduce the symmetry of the crystal 



in such a way that the positions of the Pt atoms are no 
longer completely fixed by the symmetry. The strain- 
induced forces drive them into energetically more favor- 
able positions. However, the Si atom occupies a center 
of inversion symmetry and thus Si internal displacements 
are forbidden in all cases (i.e. the strain-induced forces 
are identically zero). Symmetry also places specific re- 
strictions on the nature of the Pt displacements. The 
symmetry of strain 6, corresponding to C44, allows Pt 
internal displacements along both the [010] and [001] di- 
rections, while the inversion operation leads to the re- 
quirement that the displacements must be equal and op- 
posite for the two Pt atoms in the primitive cell. The 
symmetry of strain 4, corresponding to cn, is the same 
as the symmetry of strain 1 and both allow internal dis- 
placements only along [001]. Once again the presence of 
inversion requires that the displacements of the two Pt 
atoms be equal and opposite. Strain 5, corresponding 
to C33, and strain 2 result in the same symmetry as the 
unstrained crystal and therefore there are no internal dis- 
placements associated with these cases, since there are no 
degrees of freedom in the internal atomic coordinates of 
the unstrained crystal. Strain 3 does lower the symmetry 
but internal displacements are still symmetry- for bidden. 
This fact, combined with the lack of displacements as- 
sociated with strain 2 means that cge is unaffected. We 
note that since strains 1 and 4 result in the same symme- 
try reduction relative to the unstrained crystal there will 
necessarily be a formal symmetry-required relationship 
between the internal displacements for these two strains. 
This relationship is obtained directly from the first prin- 
ciples calculations. 

In addition to placing restrictions on the nature of the 
internal displacements, symmetry also constrains the cor- 
responding changes in the elastic constants themselves. 
We have already seen that the values of cn and C44 are 
both allowed to change as a result of internal displace- 
ments but that C33 and cgg must both remain unchanged. 
The bulk modulus Bo is also required to be unchanged 
because it represents the crystal response to hydrostatic 
pressure, corresponding to a strain e*^^-* — ^Sij which 
preserves the full symmetry of the unstrained crystal, 
just as in the case of strains 2 and 5. The expression 
for the bulk modulus in terms of the elastic constants is 
Bo = |(2cii -I- C33 -t- 2ci2 -I- 4ci3) while the energy expres- 
sion corres pon ding to strain 2 is (cn -I- Ci2 — 4ci3 -1-2033)7^ 



(see Table VII). Our symmetry arguments have required 



that neither of these expressions can change as a result of 
internal displacements and therefore the changes in cn, 
C12, and C13 must exactly cancel from these two expres- 
sions (we have already shown in conjunction with strain 
5 that C33 cannot change) . The only way to achieve both 
cancellations is if the displacement-induced change in C13 
is identically zero and if the changes in cn and C12 are 
equal and opposite. Moreover, since cn appears as the 
sole coefficient in the energy expression corresponding to 
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strain 4 and since internal displacements can only lower 
the energy, we conclude that the value of cn must either 
decrease or remain the same. This conclusion leads to 
the seeming paradox that if cn decreases then symmetry 
requires that C12 must increase which appears to con- 
tradict the fact that internal displacements must always 
lower the energy. The resolution of this seeming paradox 
comes from the fact that it is not possible to construct 
a strain in which C12 appears as the sole coefficient in 
the expression for the strain energy. It always appears 
in conjunction with cn and we have already seen that 
cii + C12 is required by symmetry to be unchanged while 
cii — C12 can either decrease or remain unchanged. 

Our results for the six independent and non-zero elastic 
constants of a-Pt2Si are given in Table We have calcu- 
lated the elastic constants for the "frozen" configuration 
[all atoms held at the positions determined solely from 
Eq. (|^)] and with the relaxation of the strain-induced 
forces on the Pt atoms. In keeping with our general sym- 
metry arguments, we find a relaxation-induced softening 
of Cii by 4% and of C44 by 17%. In addition, C12 in- 
creases by 6% while the remaining elastic constants are 
unchanged to within numerical uncertainties. Our re- 
sults are also consistent with the symmetry requirement 
that the changes in cn and C12 be equal and opposite, 
since cn decreases by 14.8±1.5 GPa whereas C12 increases 
by 14.6±1.6 GPa. The bulk moduli calculated from the 
tetragonal elastic constants and from the fit to a Birch- 
Murnaghan equation of state are almost the same, giving 
a consistent prediction of Bq = 235 GPa. As required by 
symmetry, the bulk modulus has the same value in the 
frozen and relaxed calculations. 

The requirement that the crystal be stable against any 
homogeneous elastic deformation places restrictions on 
the elastic constants, just as in the cubic case. For tetrag- 
onal crststals these mechanical stability restrictions are as 
foUowsEfl 

(cii - C12) > 0, (cii + C33 - 2ci3) > 0, 
cii > 0, C33 > 0, C44 > 0, C66 > 0, (13) 

(2cii-|-C33-|-2ci2-f 4ci3) > 0. 

The elastic constants in Table ^ satisfy all of the con- 
ditions in Eq. (^3|). In particular, C12 is smaller than 
cii and C13 is smaller than the average of cn and C33. 
The stability conditions again lead to restrictions on the 
magnitude of Bq . We first rewrite Bq as 

-Bo = ^ [6cii -I- 3c33 - 2(cii - C12) - 2(cii + C33 - 2ci3)] • 

(14) 



for cii and C33 we obtain 



Using Eq. ( |l4| ) and the first two inequalities in Eq. (13) 
we obtain the following result 



1 . 

Bq < -(2cii + C33), 



(15) 



-Bo > 2(ci2 + 2ci3), 



(16) 



that the bulk modulus must be larger than the weighted 
average of C12 and C13. 

The stability restrictions do not tell us anything fur- 
ther about the relative magnitudes of the various elastic 
constants. For example, we find a small value of C44 in 
comparison to cee which means that the tetragonal unit 
cell is more easily deformed by a pure shear about the a- 
or b-axis in comparison to the c-axis. We also find that 
overall the elastic constants of a-Pt2Si are much closer 
to those of pure Pt than pure Si. In particular cn and 
C33 are similar in magnitude to cn in Pt, but all of these 
constants are approximately twice the value of cn in Si. 
Similarly C12 has approximately the same magnitude for 
both Q;-Pt2Si and Pt, although C13 in the silicide is about 
30% smaller. However, C12 in Si is a factor of four smaller. 
Conversely, the value of C44 is similar in magnitude for 
all three materials, with cqq in the silicide being a factor 
of two larger. The bulk modulus in the silicide is about 
20% smaller than in Pt but still more than a factor of 
two larger than in Si. The connection between the mag- 
nitudes of the various elastic constants and the chemical 
bonding has been explored in detail in a separate study.EZI 

In addition to the relaxed elastic constants we also ob- 
tained the values of the dimensionless parameters cf^'* 
which determine the magnitudes of the Pt internal dis- 
placements themselves, 



Au 



(I) 



(17) 



that the bulk modulus must be smaller than the weighted 
average of cn and C33. Similarly, by substituting instead 



where a is the lattice constant, / = 1,4,6 corresponds 
to the strains with symmetry-allowed internal displace- 
ments, and i is the Cartesian index (i = 2,3 for J = 6 
and i = 3 for 7=1,4). In principle there could be con- 
tributions to the Am^^-* which are of higher order in 7, 
but since we are only considering the second order elas- 
tic constants the strain energy is only expanded to second 
order in 7 [see Eqs. (||) and (||)], or equivalently, to second 
order in the total displacements, Ui + Aui. Thus we need 
only consider the linear term in Eq. (lO) in the present 

context. The calculated values of the Q parameters are 
listed in Table XI along with the displacement param- 
eter associated with C44 in Si. We note that although 
displacements are allowed along both the [010] and [001] 
directions in the case of the a-Pt2Si strain 6 (correspond- 
ing to C44), the displacenients along [001] are found to be 
zero to linear order in 70 As indicated above, we expect 
the internal displacements for strains 1 and 4 to exhibit a 
symmetry-required relationship and in keeping with this 
expectation we find that Cs^"^' is almost exactly three 
times larger than C3^~^\ the difference likely being due 
to numerical uncertainty. 
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E. PtSi 



foUowsB 



Nine independent strains are necessary to compute 
the elastic constants of orthorhombic PtSi. We first 
performed calculations of the elastic constants with 
the internal structural parameters Mpt/Si f^nd wpt/Si 
held "frozen" at their self-consistent equilibrium va lues . 
These results are listed in the second column of Table XI] . 



The E{'^) curves are well-fitted by third-order polynomi- 
als in 7, as can be seen from the small standard errors 
in the calculated Cy . The value of Bq obtained from the 
elastic constants, Bq = |(cii -I- C22 + C33 -I- 2ci2 -I- 2ci3 + 
2C23), agrees reasonably well with the one which was de- 
termined from the calculation of the lattice constants. 
This is not surprising since the lattice constant calcula- 
tions were also performed with frozen atomic degrees of 
freedom. 

As expected, the equilibrium atomic positions are not 
independent of the shape and size of the unit cell — 
similarly to the case of C44 in Si as well as cn, C12, and 
C44 in Q;-Pt2Si, as discussed above. The relaxed elas- 
tic co nsta nts for PtSi are listed in the last column of 
Table XI]. For example, C44 drops from 141.3 GPa to 
100.1 GPa when all of the atoms are relaxed. We find 
that the Si atoms adjust to this shear about the a-axis 
by moving mainly along the b-axis. Most of the other 
elastic constants decrease by 10-20%, except C23 which 
increases by 8% and C12 which remains approximately 
unchanged. In analogy with the case of C12 for tetrago- 
nal a-Pt2Si, we note that for an arbitrary strain C12, C13, 
and C23 in orthorhombic PtSi never appear isolated but 
always occur in combination with other elastic constants 
in the expression for the second-order change in the total 
energy [Eq. (^)]. These particular elastic constants are 
therefore not required to decrease when relaxation is in- 
cluded, even when the energy is lowered. Conversely, the 
remaining six elastic constants are required to decrease 
when relaxation lowers the energy because strains can 
be constructed for which each appears as the isolated 
coefficient of the only contribution to the second-order 
change in the energy. In the case of PtSi, the additional 
relaxation of the internal degrees of freedom leads to a 
significant softening of the elastic constants which must 
also be taken into account in determining the bulk mod- 
ulus. Therefore, we predict the bulk modulus of PtSi to 
be 198 GPa, which is 6% lower than the value of 210 GPa 
determined in our frozen-configuration total-energy min- 
imization. Although we have calculated the changes in 
the elastic constants when the internal atomic degrees of 
freedom are allowed to relax, we have not explicitly ex- 
tracted the corresponding internal displacement param- 
eters as we did for pure Si and a-Pt2Si. 

Mechanical stability leads to restrictions on the elas- 
tic constants, which for orthorhombic crystals are as 



(cii + C22 - 2ci2) > 0, (cii + C33 - 2ci3) > 0, 
(c22 + C33 - 2C23) > 0, 
cii > 0, C22 > 0, C33 > 0, (18) 

C44 > 0, C55 > 0, C66 > 0, 

(cii + C22 + C33 + 2ci2 + 2ci3 + 2C23) > 0. 



The elastic constants in Table XII satisfy all of these con- 
ditions and in particular, C12 is smaller than the average 
of cii and C22, C13 is smaller than the average of cn and 
C33, and C23 is smaller than the average of C22 and C33. 
As in the case of Q!-Pt2Si, we can obtain restrictions on 
the magnitude of Bq, 

i(ci2 + Ci3 -I- C23) < Bo < ^(Cii + C22 + C33), (19) 

that the bulk modulus must be smaller than the average 
of cii, C22, and C33 but larger than the average of C12, 
ci3, and C23. 

We again find that overall the elastic constants of PtSi 
are much closer to those of pure Pt and a-Pt2Si as com- 
pared to pure Si. In detail, we find that cn, C22, and C33 
are approximately 10% smaller on average in PtSi than 
in a-Pt2Si, and that Bq is 16% smaller. In addition, C12, 
C13, and C23 for PtSi are close in magnitude to C13 for 
ck-Pt2Si, which we saw was about 30% smaller than C12 
in both pure Pt and a-Pt2Si. Finally, C44, C55, and cge 
for PtSi are similar to the values of C44 in all three of the 
other materials, but still approximately a factor of two 
smaller than cge in a-Pt2Si. 



F. Trends in the elastic constants 

The trends of the elastic constants as a function of 
the atomic percent Pt in all four materials are plotted 
in Fig. ^. Each of the curves corresponds to an average 
of a different class of elastic constants, while the sym- 
bols show the values of the individual elastic constants 
themselves. As we saw in Eqs. ([T5|), (p^), and (p^), me- 
chanical stability requires that Bq be larger than the av- 
erage of Cii, C22, and C33 but smaller than the average 
of C12, Ci3, and C23 (note that in the case of a-Pt2Si the 
appropriate averages are ■|(2cii + C33) and -^(012 -I- 2ci3) 
because cn = C22 and C13 = C23 for tetragonal crystals). 
This stability requirement is reflected in the top three 
curves in Fig. ||. We also see that these three curves 
each increase monotonically from Si to Pt and we note 
that all three classes of elastic constants represented by 
these curves correspond to strains in which the volume 
is not fixed. Conversely, the two lower curves labeled 
(cii — ci2)/2 and C44 correspond to the two classes of 
elastic constants in which the strains are strictly volume- 
conserving [in the case of PtSi the lowest solid-line curve 
and large open circles correspond to elastic constant com- 
binations j(cii + C22 - 2ci2), |(cii + C33 - 2ci3), and 



7 



j{c22 + C33 — 2C23)]. We see that in this case the two 
sets of averages are approximately constant as a function 
of atomic percent Pt. The significance of this difference 
in the trends of volume-conserving versus non-volume- 
conserving elastic constants is connected to the nature 
of the chemical bonding in these materials and has been 
addressed in a separate study.E3 



IV. ELECTRONIC STRUCTURE 



In a typical good metal such as Pt, the optical absorp- 
tion at low energies is dominated by the free-electron- 
like Drude contribution. However, we saw from Fig. ^ 
that in the case of the silicides they have a low DOS at 
the Fermi level and consequently are poor metals. In 
this circumstance the Drude contribution will be greatly 
reduced and therefore the presence of many low-energy 
splittings in the bands near the Fermi level may result 
in an inter-band contribution to the optical absorption 
which is significant even at low energies in the infrared 
range. 



The self-consistent calculations for the spin-orbit-split 
energy bands of a-Pt2Si andpPtSi were performed using 
the WIEN97 implementationEp-|(^ the linear augmented 
plane wave (LAPW) methodBj-ES The local density ap- 
proximation was used with the-iexchange-correlation po- 
tential of Perdew and Wang.cll The effects of the spin- 
orbit interaction .wrgre included in a second-order varia- 
tional procedureB^H In the self-consistency cycles ap- 
proximately 120 irreducible k-points (1000 k-points in 
the full BZ) wepe, used in the modified tetrahedron 
method of BlochLcJ The energy cutoff used for the plane- 
wave expansion was kmax = 4.16 a.u. resulting in a well 
converged basis set of about 105 basis functions per atom. 
The experimental values of the lattice constants and in- 
ternal structural parameters from Tables and |T| were 
used in all cases. For the purpose of calculating the 
density of states (DOS) we again used the tetrahedron 
method but with unshifted k-point meshes which in- 
cluded the r point. In the case of fee Pt, cubic-diamond- 
phase Si, and a;-Pt2Si we used a 32x32x32 mesh cor- 
responding to 897 irreducible k-points for the two cubic 
materials and 2393 for the silicide. For PtSi we used 
a 16x24x16 mesh, yielding 1053 irreducible k-points 
points. 

The total DOS for all four materials is shown on the 
same scale in Fig. ^ Although both silicides are metals 
with a non-zero DOS at the Fermi level, they are found 
to be poor metals since the DOS is much smaller than 
in the case of pure Pt which is a good metal. In pure 
Pt the Fermi level lies near the top but still within the 
large DOS features of the d— band but in both sihcides 
the Fermi level lies above these large peaks. In addi- 
tion, both silicides exhibit a peak in the DOS at around 
— 10 eV which arises from the Si s-orbitals. We note that 
the basic features of the electronic structure, as reflected 
in the total DOS, do not appear to differ very much be- 
tween the two silicides. The origin of the various features 
in the-PtSi DOS has been discussed in detail by Franco 
et alB 

We have also calculated the spin-orbit-split energy 
bands near the Fermi level for the two silicides, as shown 
in Fig. H In both cases we see that there are a size- 
able number of small-energy splittings between different 
bands throughout the full BZ. These small splittings are 
of direct interest with regard to low-energy inter-band 
transitions which contribute to the optical absorption. 



V. SUMMARY 



We have carried out an extensive first principles study 
of two room-temperature stable Pt silicides, tetragonal 
a-Pt2Si and orthorhombic PtSi. We have determined 
the theoretical equilibrium structural parameters and co- 
hesive energies for both silicides, as well as pure fee Pt 
and pure cubic-diamond-phase Si. In particular, we have 
carried out a large number of calculations in order to 
minimize the total energy with respect to the two lat- 
tice constants in tetragonal a-Pt2Si and the three lattice 
constants and four internal structural parameters of or- 
thorhombic PtSi. Our calculated structural parameters 
for all four materials are in good agreement with experi- 
mental data, validating the method we have used. 

A major portion of our effort here has been directed 
at the elastic constants in the two silicides. All of the 
independent, non-zero elastic constants (6 for a-Pt2Si, 
9 for PtSi, and 3 each for the two cubic materials) have 
been calculated from first principles. The silicide calcula- 
tions required extensive relaxation of the internal degrees 
of freedom, especially in the case of the low symmetry 
structure of orthorhombic PtSi. Comparing the elastic 
constants obtained with and without relaxation we find 
that relaxation induces significant changes in the mag- 
nitudes of many of the elastic constants. In addition, 
we have explicitly determined the dimcnsionless internal 
displacement parameters for the three strains in a-Pt2Si 
for which they are non-zero. We also note that the value 
of C44 in pure fee Pt was found to be extremely sensitive 
to the number of k-points, much more so than any of the 
other elastic constants we calculated. This sensitivity 
results from the large number of van Hove singularities 
close to the Fermi level. 

We have investigated the trends in the calculated elas- 
tic constants, both the trends within a given material 
as well as between materials. The requirement of me- 
chanical stability places specific restrictions on the rela- 
tive magnitudes of some of the elastic constants within a 
given material, including, for example, a restriction on 
the bulk modulus that Bq < |(cii + C22 + C33) and 
Bq > |(ci2 -f- Ci3 -I- C23). With regard to the trends 
among the four materials, we find that in the metals the 
elastic constant expressions which correspond to volume- 
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conserving strains are always smaller than those which 
correspond to strains which do not conserve volume. This 
also turns out to be true in Si with the exception of 
ci2 which is less than C44. However, the difference in 
magnitudes between volume-conserving and non- volume- 
conserving elastic constants is largest on average in Pt 
and gets smaller in the progression Pt ck-Pt2Si PtSi 
—>■ Si. In general, the volume-conserving elastic constants 
have similar magnitudes in all four materials while the 
non- volume-conserving elastic constants follow this same 
progression. In particular, the bulk modulus is found to 
be a very nearly linear function of the atomic percentage 
of Pt. Klepeis et a?£3 have studied the close connection 
between the various trends in the elastic constants and 
the chemical bonding in the Pt silicides. 

The calculated electronic structure demonstrates that 
the two silicides are poor metals with a low density of 
states at the Fermi level, and consequently we expect that 
the Drude component of the optical absorption should be 
much smaller than in good metals such as pure Pt. In 
addition, we find a large number of small-energy differ- 
ences between various bands near the Fermi level in the 
calculated spin-orbit-split band structure for the two sili- 
cides. These two circumstances suggest that it may be 
important to include the interband contribution to the 
optical absorption as well, even in the infrared region. 
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FIG. 1. Conventional unit cells of (a) body-centered tetrag- 
onal a-Pt2Si and (b) orthorhombic PtSi. The relevant lattice 
constant distances are illustrated in both cases. 
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FIG. 2. Trends in the elastic constants as a function 
of atomic percent Pt for pure cubic-diamond-phase Si, or- 
thorhombic PtSi, tetragonal Q-Pt2Si, and fee Pt. The differ- 
ent curves correspond to the average values of different classes 
of the individual elastic constants, as specified in the legend. 
For example, in the case of the dotted-line curve labeled as 
ci2, the line passes through ^(ci2 + Ci3 + C23) in the case of 
PtSi and through |(ci2 + 2ci3) for Q-Pt2Si (ci3 = C23 for 
tetragonal crystals), while the open squares show the actual 
values of C12, C13 and C23, as appropriate for each material. 
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FIG. 3. Total density of states plotted on the same scale 
for the four materials considered in this work, ordered by 
decreasing Pt content from top to bottom. Only Si is a semi- 
conductor whereas Pt and the two silicides are metals. The 
valence band maximum in Si is labeled as the Fermi level Ep. 
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FIG. 4. Spin-orbit split energy bands near the Fermi level 

for (a) a-Pt2Si and (b) PtSi. The primary feature of inter- 
est in both plots is the relatively large number of low-energy 
splittings between various bands near the Fermi level. 
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TABLE I. Structural characterization of the materials studied in this work. The column labeled "Structure" gives the 
Strukturbericht designations. The columns "Space Group" (name and number), "Site" (multiplicity and Wyckoff letter), and 
"1st position" follow Ref. Wq. See Table for the values of the internal structural parameters u and v. 



Material 


Structure 


Space group 




Site 




1st position 


Ref. 


Pt 


Al 


Fm3m 


225 


Pt 


4a 


0,0,0 


P 


a-PtaSi 


L'2t 


lA/mram 


139 


Pt 


Ad 


i i 

^' 2 ' 4 


§ 










Si 


2a 


0,0,0 




PtSi 


B31 


Pnma 


62 


Pt 


4c 


upt, J,VPt 












Si 


4c 


USi, l,VSi 




Si 


A4: 


FdSm 


227 


Si 


8a 


111 

8 ' 8 ' 8 






TABLE IL Theoretical and experimental lattice constants (in a.u.) and bulk moduli Bo (in GPa). The theoretical Bo were 
determined from a fit to a Murnaghan equation of state for Pt and Si, and to a four-term Birch-Murnaghan equation of state 
for a-Pt^Si and PtSi. Bo listed here for PtSi was calculated for fixed values of the internal structural parameters (see also 
Table XII). The experimental value of the bulk modulus for Pt is extrapolated to K, whereas all of the other experimental 
numbers are given for room temperature. 



Material 




ao 


bo 


Co 


Bo 


Ref. 


Pt 


thcor. 


7.403 






287.8 






exp. 


7.415 






288.4 


a-Pt2Si 


theor. 


7.407 




11.241 


233.5 




exp. 


7.461 




11.268 




i 


PtSi 


theor. 


10.583 


6.774 


11.195 


210.0 






exp. 


10.539 


6.778 


11.180 




Si 


theor. 


10.22 






95.9 




exp. 


10.26 






98.8 





TABLE III. Calculated and experimental internal structural parameters of PtSi. The parameters u and v are the same 
as those specified in Table |. These parameters were calculated both for the experimental lattice constants as well as the 
self-consistent theoretical lattice constants (Table |l|). The experimental internal parameters for PtSi are also listed. 



PtSi 


upt 


vpt 




"Si 


At expt. lattice constants 


0.9981 


0.1915 


0.1777 


0.5845 


At theor. lattice constants 


0.9977 


0.1919 


0.1782 


0.5841 


Experimentlij 


0.9956 


0.1922 


0.177 


0.583 



13 



TABLE IV. Parameters describing the basis used in the FPLMTO calculations. R^t is the muffin-tin radius in a.u., Lmax 
is the upper limit on the angular momentum expansion of the smoothed Hankel functions about a given atomic site, /("max is 
the order of the biorthogonal polynomials used in this expansion, iiam is the smoothing radius in a.u., and —K^ is the decay 
energy. The total number of basis functions per atom is 17 for Ft and 13 for Si. See Ref. |l^ for a more complete description 
of the parameters. 



Atom 




-^max 


-^max 


L -DiOCK 


Basis 

-n,sm 


,J2 


Pt 


2.2 


3 


4 


s 


2.865 


-1.02 










p 


2.130 


-1.08 










p 


1.302 


-1.53 










d 


1.000 


-0.89 










d 


2.123 


-0.42 


Si 


2.1 


3 


5 


s 


1.908 


-1.296 










p 


1.627 


-0.302 










d 


1.601 


-1.496 










s,p 


2.200 


-2.000 



TABLE V. Cohesive energies iJcoh and heats of formation A//f in eV/atom. The experimental standard heat of formation is 
given for T = 298.15 K whereas theoretical values are valid for K and do not contain any corrections for zero-point vibrations. 
The experimental cohesive energies for Pt and Si also correspond to K. 



Material 




-Bcoh 




Ref. 


Pt 


thcor. 


7.27 









exp. 


5.84 





§ 


a-PtaSi 


theor. 


7.24 


-0.65 






exp. 


6.08 


-0.64 


PtSi 


theor. 


6.93 


-0.67 






exp. 


5.85 


-0.62 


Si 


theor. 


5.23 





§ 




exp. 


4.63 






TABLE VL Parameterizations of the three strains used to calculate the three elastic constants of cubic Pt and Si (also used 
in Ref. ^). The energy expressions were obtained from Eq. (^). Strains 7 = 1 and / = 3 are strictly volume-conserving to 
all orders in the strain parameter 7. If we restrict ourselves to linear order only then eg^' — —27 and 63'^'' — 0, with volume 
conservation preserved to linear order as well. 



Strain / 


Parameters (unlisted = 0) 


I\EIV to 0(7^) 


1 


ei = 62 = 7, 63 = (1 + 7)"^ - 1 


3(cii - Ci2)7^ 


2 


61 = 62 = 63 = 7 


|(cii +2ci2)7^ 


3 


66 =7, 63 =7^(4-7^)"^ 


|C447^ 



14 



TABLE VII. Parameterizations of the six strains used to calculate the six elastic constants of tetragonal a-Pt2Si (taken from 

Ref. |2^) . The energy expressions were obtained from Eq. (^) . 

Strain / Parameters (unlisted = 0) /^E/V to 0(7^) 

1 ei = 27, 62 = 63 = -7 |(5cii - 4ci2 - 2ci3 + 033)7^ 

2 ei = 62 = -7, 63 = 27 (cii + C12 - 4ci3 + 2633)7^ 

3 61 = 62 = 7, 63 = -27, 66 = 27 (Cii + C12 - 4ci3 + 2633 + 2666)7^ 

4 61 = 7 icii7^ 



5 63 = 7 56337^ 

6 64 = 27 2c447^ 



TABLE VIII. Parameterizations of the nine strains used to calculate the nine elastic constants of orthorhombic PtSi. The 
energy expressions were obtained from Eq. (^). 

Strain / Parameters (unlisted = 0) /^E/V to 0(7^) 

1 61 = 7 ^cii7^ 

2 62 = 7 56227^ 



3 63 = 7 ^6337^ 



4 61 = 27, 62 = -7, 63 = -7 1(4611 - 4ei2 - 4ei3 + 622 + 2e23 + 633)7^ 

5 61 = -7, 62 = 27, 63 = -7 I (611 - 4612 + 2ei3 + 4622 " 4623 + 633)7^ 

6 61 = -7, 62 = -7, 63 = 27 i(6ii + 2612 - 4613 + 622 - 4623 + 4633)7^ 
I? 12 

7 64 = 7 

8 65 = 7 ^6557^ 

9 66 = 7 ^6667^ 



TABLE IX. Elastic constants of Pt and Si. Calculations were carried out at the theoretical self-consistent lattice constants 
of apt = 7.403 a.u. and asi = 10.22 a.u. The theoretical value of 644 in parentheses for Si is the "frozen" value obtained 
without allowing for internal relaxation. The bulk modulus is calculated from the elastic constants as Bo ~ ^(cii + 2ci2). In 
parentheses we give Bo from the fit to a Murnaghan equation of state. Experimental values are extrapolated to K. All values 
are in units of GPa. 





Pt Theory 


Pt Expt.S 


Si Theory 


Si ExptH 


Cll 


346.8±0.5 


358 


163.45±0.03 


165 


612 


262.7±0.3 


254 


62.13±0.02 


63 


644 


87.5±0.3 


77 


79.85±0.02 (108.6) 


79.1 


Bo 


290.8±0.3 (287.8) 


288.4 


95.90±0.02 (95.9) 


97.0 
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TABLE X. Elastic constants of Q-Pt2Si. Calculations were performed at the theoretical self-consistent lattice constants 
(Table 0). "Frozen" refers to fixed atomic positions, whereas "relaxed" indicates that a relaxation of the atomic positions 
was carried out. Parentheses denote values where no internal relaxation was necessary because of symmetry constraints (small 
variations in these values come from using a slightly more stringent convergence criterion on the energy). The bulk modulus is 
calculated from the elastic constants as Bo — ^(2cii +C33 +2ci2 +4ci3). B^^'^'^^ is from a Birch-Murnaghan fit. No experimental 
data is available. All values are in units of GPa. 

Qf-Pt2Si frozen relaxed 

cTi 347.2±1.2 332.4±0.9 

C33 297.5±0.5 (298.0±0.4) 

ci2 225.0±1.2 239.6±1.0 

ci3 169.3±0.9 (169.4±0.8) 

C44 75.4±0.3 62.7±0.5 

C66 169.5±5.2 (169.3±5.2) 

Bo 235.4±0.6 (235.5±0.5) 

B^"'''^ 233.5 — 



TABLE XL Dimensionless internal displacement parameters for Si [Eq. (p^, strain 3 from Table ^| and a-Pt2Si [Eq. (p^, 
strains 1, 4, and 6 from Table VII |. In the case of strain 6 for Q-Pt2Si displacements were calculated a ong both the [010] and 
[001] directions, whereas displacements only along the [001] direction are allowed by symmetry for strains 1 and 4. In addition, 
there is a strict symmetry-required relationship between the displacements in strains 1 and 4 (see text). Experimental data is 
available only for Si. 



Material 



Strain I 



Theory 



Experiment 



Ref. 



Si 

a-Pt2Si 



3 
1 

4 

6 [010] 
6 [001] 



0.53 
0.22 
0.074 
-0.12 
0.00 



0.54 



TABLE XII. Elastic constants of PtSi calculated at the theoretical self-consistent lattice constants (Table ||). The second 
column shows the elastic constants obtained when the internal structural parameters were held fixed at their theoretical 
self-consistent values (Table In the third column these were allowed to relax. The bulk modulus is calculated from the 

elastic constants as Bo = |(cii -I- C22 + C33 -I- 2ci2 -I- 2ci3 -I- 2C23). Bq"^''^ is from a Birch-Murnaghan fit obtained with frozen 
values of the internal structural parameters. No experimental data is available. All values are in units of GPa. 



PtSi 


frozen 


relaxed 


Cll 


327.5±1.2 


298.2±1.2 


C22 


313.8±0.0 


269.3±0.8 


C33 


345.9±0.1 


308.0±0.6 


C12 


157.7±0.6 


156.4±0.8 


Cl3 


162.9±0.6 


132.2±0.7 


C23 


153.4±0.1 


165.1±0.6 


C44 


141.3±0.3 


100.1±0.4 


C55 


113.1±0.1 


104.5±0.1 


C66 


74.2±0.2 


66.3±0.4 


Bo 


215.0±0.2 


198.1±0.3 


D Birch 

-fc>o 


210.0 
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